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Abstract 

We propose a general integrable lattice system involving some free 
parameters, which contains known integrable lattice systems such as 
the Ablowitz-Ladik discretization of the nonlinear Schrodinger (NLS) 
equation as special cases. With a suitable choice of the parameters, it 
provides a new integrable space-discretization of the derivative NLS 
equation known as the Chen-Lee-Liu equation. Analogously to the 
continuous case, the space-discrete Chen-Lee-Liu system possesses a 
Lax pair and admits a complex conjugation reduction between the two 
dependent variables. Thus, we obtain a proper space-discretization of 
the Chen-Lee-Liu equation defined on the three lattice sites n — 1, 
n, n + 1 for the first time. Considering a negative flow of the dis¬ 
crete Chen-Lee-Liu hierarchy, we obtain a proper discretization of 
the massive Thirring model in light-cone coordinates. Multicompo¬ 
nent generalizations of the obtained discrete equations are straight¬ 
forward because the performed computations are valid for the general 
case where the dependent variables are vector- or matrix-valued. 
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1 Introduction 


There exists no systematic method of obtaining a proper discretization for a 
given integrable partial differential equation (PDE). For some scalar PDEs 
such as the sine-Gordon equation u xt = sin u M, elementary auto-Backlund 
transformations and the associated nonlinear superposition formula based 
on Bianchi’s permutability theorem can provide proper discretizations of the 
original continuous equations. However, this idea does not apply directly to 
complex-valued PDEs involving the operation of complex conjugation, such 
as the nonlinear Schrodinger (NLS) equation i q t + q xx ± 2 q 2 q* = 0, which can 
be naturally obtained from two coupled PDEs by imposing a complex conju¬ 
gation reduction on the two dependent variables [3]. This is because elemen¬ 
tary auto-Backlund transformations do not, in general, maintain the complex 
conjugation reduction [SHE]- One can consider more elaborate auto-Backlund 
transformations by composing two or more elementary auto-Backlund trans¬ 
formations so that the net result can maintain the complex conjugation 
reduction between the two dependent variables. However, such composite 
auto-Backlund transformations generally involve an indefinite integral [Q] or 
a square-root function with an indefinite sign [TOUTT] (also see (7.18) in [12] 
or (4.17) in [L3j), so they do not directly provide proper discretizations of 
the original PDEs, which can be written in local form and define the time 
evolution uniquely in the discrete setting. 

Because proper discretizations of integrable PDEs involving the operation 
of complex conjugation cannot be derived in a systematic manner, it is more 
productive and instructive to obtain such discretizations on a case-by-case 
consideration. In this paper, we propose a new proper space-discretization of 
an integrable derivative NLS equation i q t + q xx + i qq*q x = 0, usually referred 
to as the Chen-Lee-Liu equation [14| (also see [IS]). Actually, we already 
obtained an integrable space-discretization of the Chen-Lee-Liu equation in 
our previous paper [16], but it depends on five lattice sites and is rather 
complicated. For the Chen-Lee-Liu system in the nonreduced form, i.e., two 
coupled PDEs for q and q* wherein q and q* are independent functions not 
related by a complex conjugation, a fully discrete analog was derived in [8] 
(see also the continuous-time flows in mmy, a correspondence between a 
fully discrete system and continuous-time flows is briefly described in EH- 

We can obtain a proper space-discrete Chen-Lee-Liu equation defined 
on three lattice sites n — 1, n, n + 1 by constructing its Lax-pair representa¬ 
tion. As the spatial part of the Lax-pair representation, we consider a new 
discrete spectral problem implied by a binary Backlund-Darboux transfor¬ 
mation for the continuous Chen-Lee-Liu equation or, equivalently, that for 
the Ablowitz-Ladik lattice (an integrable discrete NLS equation [221123 ] ) pro- 
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posed in [24,25J- Note that the Chen-Lee-Liu hierarchy and the Ablowitz- 
Ladik hierarchy are different facets of the same object and in a sense equiv¬ 
alent j26jL22j (also see [28([29]). In contrast to the usual formulation of the 
Backlund-Darboux transformations [5HZU30H32] , we do not express two un¬ 
known functions appearing in the Darboux matrix explicitly in terms of linear 
eigenfunctions of the original spectral problem. Rather, we consider them 
as new dependent variables in the discrete spectral problem, which can be 
related to each other through a complex conjugation reduction. Then, we 
associate the discrete spectral problem with a suitable time-evolution equa¬ 
tion, which comprises the Lax-pair representation; the compatibility condi¬ 
tion provides a rather general lattice system involving some arbitrary pa¬ 
rameters. Depending on the choice of the parameters, it provides a proper 
space-discretization of the Chen-Lee-Liu equation, as well as other integrable 
lattice systems such as the Ablowitz-Ladik lattice [22423] . 

This paper is organized as follows. In section 2, we introduce a discrete 
spectral problem and associate it with a suitable isospectral time-evolution 
equation. Then, the compatibility condition for this Lax-pair representation 
provides a new integrable lattice system; with a suitable choice of the param¬ 
eters therein, we obtain a proper space-discretization of the Chen-Lee-Liu 
equation. In section 3, we change the time part of the Lax-pair represen¬ 
tation to obtain a proper discretization of the massive Thirring model in 
light-cone (or characteristic) coordinates [TI31 [TTb] . Section 4 is devoted to 
concluding remarks. Throughout the paper, we perform the computations 
for the most general case where the dependent variables are (rectangular) 
matrix-valued, so multicomponent generalizations of the obtained equations 
are straightforward. 

2 Space-discrete Chen—Lee—Liu equation 

In this section, we propose a new integrable lattice system, which contains a 
proper space-discretization of the Chen-Lee-Liu equation as a special case, 
through a Lax-pair representation. In subsection 12.11 we introduce a new 
discrete spectral problem, which gives the spatial part of the Lax pair. In 
subsection E2 we associate it with an isospectral time-evolution equation, 
which is the temporal part of the Lax pair. 
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2.1 Backlund—Darboux transformation as a discrete 
spectral problem 

The Ablowitz-Ladik spectral problem [221123] can be generalized to a block- 
matrix form [36] (also see [TBlfST] and references therein) as given by 

" ^ 1,771+1 

V I / 2,m+l 

Here, ( is a constant spectral parameter and Q m and R m are, respectively, 
li x I 2 and I 2 x l\ (generally rectangular) matrices. Thus, the square matrix 
above is partitioned as an ( l\ + 12 ) x (l 1 + 12 ) block matrix; for notational 
brevity, we omit the index of each unit matrix I to indicate its size. 

Among an infinite set of isospectral time-evolution equations compati¬ 
ble with the Ablowitz-Ladik spectral problem (12.111 (cf. [SBiES]), the most 
fundamental ones are 
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(2.3) 


Here, the subscript tj denotes time differentiation and the symbol italic O 
represents a zero matrix. The corresponding equations of motion are derived 
from the compatibility condition for the overdetermined linear system, (12.ip 
and (12.21) or (12.31) . i.e., 


Qm,ti Qm+ 1 T Qm+lRmQm O, 
Rm,ti T Rm— 1 RmQmRm— 1 O, 


(2.4) 


and 


Qm,t —i d" Qm —1 QmRmQm —1 Oi 
Rm,t— i Rm-\- 1 H“ O. 


(2.5) 


These are the two elementary flows of the Ablowitz-Ladik hierarchy and 
a suitable linear combination of them together with the trivial zeroth flow 
provides the Ablowitz-Ladik discretization of the NLS system [221 HU] . Note 
that the use of O , instead of 0, on the right-hand side of the equations implies 
that the dependent variables are matrix-valued. 
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A crucial fact noticed by Vekslerchik [2BJI27J (also see Barashenkov- 
Getmanov [2S1S2]) is th a t th e Ablowitz-Ladik hierarchy is in a sense equiv¬ 
alent to a derivative NLS hierarchy called the Chen-Lee-Liu hierarchy US- 
This is clear from the observation that the isospectral equation (12.2|) with 
ti —> x (or (12.3)) with t_i —> x ) has essentially the same form as the spec¬ 
tral problem associated with the continuous Chen-Lee-Liu hierarchy [4T1 EFT] 
(see gHllS] for the matrix case); note that this spectral problem can be 
rewritten in a traceless form using a simple gauge transformation [28]. Ac¬ 
tually, the Chen-Lee-Liu spectral problem was first presented for the mas¬ 
sive Thirring model in light-cone (or characteristic) coordinates [33H35] (also 
see (46.,;4?j for the case of laboratory coordinates), which is the first nega¬ 
tive flow of the Chen-Lee-Liu hierarchy [44Jj48] ! t-he compatibility condition 
for the overdetermined linear system, (12.2p and ()2.3j) for any fixed value of 
77i, indeed provides the equations of motion for the massive Thirring model 
((E3D for Q m _ Ml and R mM and (I2A)) for Q m>t _ 1 and R m - l.t.J. Conversely, 
we can regard the Ablowitz-Ladik spectral problem (12.11) as a Backlund- 
Darboux transformation for the continuous Chen-Lee-Liu hierarchy [2411ZT31 ; 
then the elementary flow (12.41) (or (12.51) 1 of the Ablowitz-Ladik hierarchy 
can be considered as an auto-Backlund transformation for the continuous 
Chen-Lee-Liu hierarchy. These items of information imply that a suitably 
chosen Backlund-Darboux transformation for the Ablowitz-Ladik hierarchy 
can provide a proper discretization of the Chen-Lee-Liu spectral problem. 

In fact, Backlund-Darboux transformations for the Ablowitz-Ladik hier¬ 
archy have already been considered in a number of papers (see, e.g., [2411321 
091 - 50 ]); those results are closely related to integrablc time-discretizations of 
the flows of the Ablowitz-Ladik hierarchy [2511271l38ll40ll5Tll52] . Noting that 
the Ablowitz-Ladik hierarchy is invariant under the rescaling Q m —>- kQ m , 
R rn —y k~ l R m , we can introduce a slightly generalized version of the binary 
Backlund-Darboux transformation in the form (cf. [241125) 153]): 
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^2 ,m 


Here, a, (3, 7 and S are arbitrary parameters and u m and v m are intermedi¬ 
ate potentials that can be expressed in terms of linear eigenfunctions of the 
Ablowitz-Ladik spectral problem (12. ip at special values of the spectral pa¬ 
rameter <C; however, this information is not necessary for the purpose of this 
paper and we will treat u m and v m as new independent unknowns. Alter¬ 
natively, we can consider (12. 6 [) as a generalized binary Backlund-Darboux 
transformation for the isospectral problem ([2.2[) or (12.31) associated with 
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the Chen-Lee-Liu hierarchy. The binary Backlund-Darboux transformation 
for the Chen-Lee-Liu hierarchy (or, more specifically, the massive Thirring 
model) was apparently first considered by P. I. ffolod (also written as Golod, 
preprint in Russian, Kiev, 1978) in a somewhat preliminary form; the result 
can be found in more accessible papers [Ml 156] (see [28ll29ll5Tll58] for further 
developments). 

Because the Ablowitz-Ladik hierarchy is invariant under a space trans¬ 
lation m —> m + l, Backlund-Darboux transformations can always be com¬ 
bined with a space translation; for instance, using (12.6[) and (12.1 [) . we can 
compute 

d(i,m+i 
^2,m+1 


or 


1 


and consider it as a new Backlund-Darboux transformation, which has the 
same ^-dependence as that usually employed in the time-discretizations of the 
Ablowitz-Ladik flows • Note also that an overall factor of (I2.6|l 

is nonessential, so we can multiply the right-hand side of (12.6|) by ( or l/£ to 
obtain more familiar ^-dependence that appears in the time-discretizations 
of the Ablowitz-Ladik flows. 

Now, we freeze the lattice index m for the Ablowitz-Ladik hierarchy and 
reinterpret the Backlund-Darboux transformation (12.6[) as defining a new 
discrete spectral problem: 


^l.n+l 

V I / 2,n+l 

where n can be intuitively understood as the number of iterations of the 
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Backlund-Darboux transformation. Here, the Lax matrix L n is given by 
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The parameters a, /3, 7 and d can be varied at each application of the 
Backlund-Darboux transformation, so they can be arbitrary functions of the 
discrete independent variable n; however, for simplicity we consider them 
as constants. Each of the four equivalent expressions in (12. 8 p has its own 
advantages. 


2.2 Isospectral time-evolution equation 


To compose a Lax pair, we associate (12.71) with a suitable isospectral time- 
evolution equation, 


'hpn 

^ 2 ,n 


= M n 


*l,n 

^ 2 ,n 


(2.9) 


The compatibility condition for the overdetermined linear system, (12.71) and 
(12.91) . is given by (a space-discrete version of) the zero-curvature equation 


L n ,t + L n M n — M n+ 1 L n — O, (2-10) 

where L n and M n comprise the Lax pair. For the Lax matrix L n in (12.81) . 
the first expression (j2.8a.[) allows us to compute L nt in such a way that the 
time-evolution equations for u n and v n can be obtained from (I2.10|) explicitly; 
note that (A"~ 1 ) t = — X~ x X t X~ l for a square matrix X. 

To obtain a proper space-discrete analog of the Chen-Lee-Liu system, we 
consider a temporal Lax matrix M n that goes well with the factorized form 
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(I2.8cp or f!2.8cip of the spatial Lax matrix L n . Thus, a natural ansatz for M n 

is 
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where F n and G n are square matrices, and c and d are arbitrary constants 
introduced for adding the trivial zeroth flow. The condition a/3 — id ^ 0 is 
assumed for the Lax pair to be nontrivial. Substituting (j2.8j) and ( 12 . lip into 
( 12 . 101 ) . we obtain recursion relations for determining F n and G n ; they are 
satisfied by setting 

F n = a (Pi I - u n _iu n ) _1 , G n = b(a5I - u n _iM n ) _1 , (2.12) 


where a and b are arbitrary constants. We can also consider the more general 
case where a and b , as well as c and d in (12. lip , are arbitrary functions of 
time t, but we do not discuss it in this paper. Then, (12.101) for (12. 8p and 
(12. lip with (12. 12p provides an evolutionary lattice system: 

u n ,t ~ a l (a/31 - u n v n ) (P'yl - u n _iv n y l (/3u n - Su n _ i) 

- b5 ( 7 M n+ i - au n ) (aSI - v n u n+l y l (a/31 - v n u n ) - (c - d)u n = O, 

< 

v n ,t ~ b5 (a{31 - v n u n ) (aSI - v n -i u n ) _ (av n - 7 u n _i) 

- a 7 (5v n+ 1 - /3v n ) ((3'yl - u n v n+1 y l (a/31 - u n v n ) + (c - d)v n = O. 

(2.13) 

In the case of /3 = a*, 5 = 7 *, b = a* and d — c*, system (I2.13p admits the 
complex conjugation reduction v n = au* n with a real constant a. To consider 
a Hermitian conjugation reduction between u n and v n , we need to rewrite 
(12 .131) . By setting c — d = —a (a/3 — jS) and using the identities 

7 (a/31 - u n v n ) (f3jI - Un-iVny 1 (f3u n - 5u n - 1 ) - (a/3 - 7 ^) u n 

S ( 7 u n au n — 1 ) ( 7 u n au n —\) v n (/d'yl u n —\v n / (f3u n 5u n — 1 ) 

= (3 ( 7 U n ~ aU n -l) (Pi I - VnUn^y 1 ( 7 SI - V n U n ) , 

7 (5v n+ 1 - /3v n ) (PiI - u n v n+ i )~ 1 (api - u n v n ) + (cc/3 - iS ) v n 
= 5 (av n+ 1 - iv n ) + (Sv n+ 1 - pv n ) (Pi I - u n v n+i y l u n (av n+1 - iv n ) 

= P ( 161 - v n u n ) (PiI - v n+ iu n y l (av n+ 1 - iv n ), 
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we can rewrite (12.131) as 


" u n ,t ~ a/3 ( 7 « n - cm„_i) (£ 7 / - v n u n _ 1 ) 1 ( 7 dl - v n u n ) 

- 65 ( 7 w n+ i - au n ) (ad I - n n w„ + i ) _1 (a/3/ - n n w n ) = O, 

x (2-14) 

v re ,t - bS (a/31 - v n u n ) (ah/ - i> n _iu n ) (av n - 7 v n _i) 

- a/3 ( 75 / - n n a n ) (£ 7 / - t>„+iit n ) _1 (av n +i - 7«n) = O. 

System (12.131) (or (12.14]) ) is a rather general system involving several free 
parameters and thus encompasses simpler lattice systems as particular (or 
limiting) cases; some of them are already known. 

• By rescaling the dependent variables and parameters as u n v n = a(3u! n v ' n , 
a = a'/a and b = W//3, omitting the prime and taking the limit a, /3 —> 0 , 
(I2.13P reduces to the matrix generalization [35] (also see [151137] and 
references therein) of the Ablowitz-Ladik lattice [223 : 

j ^n,t T ad (I Zl n V n / Zl n —\ (4 DdA.) (c d)vj n O, 

\ Vn,t + 67 (I - V n Un) Vn -1 ~ a8v n+ 1 (I - U n V n ) + (c - d)v n = O. 

This is a linear combination of the two elementary flows (12.4[) and (]2.5[) 
and the trivial zeroth flow of the Ablowitz-Ladik hierarchy. 

• By rescaling the dependent variables and parameter as u n v n = /3 r )u , n v , n 
and b = b'/f3, omitting the prime and taking the limit (3 —> 0, (I2.14|) 
reduces to the already known space-discretization of the Chen-Lee- 
Liu system mm-- 


Un.t, ad ('yUn au n — 1 ) (I v n u n — 1 ) ('~yu n j r \ au n ) (a/ 'yv n u n ') O , 

a 

jj 

v n ,t -(a/ - 7 v n u n ) ( av n - 7 ^- 1 ) - ad (I - v n+1 u n y l (av n+ i - ^v n ) = O. 

a 


In the case of scalar u n and v n , this system was previously studied 
in [HH20]; in fact, it can be considered as a continuous-time analog of 
the fully discrete Chen-Lee-Liu system proposed by Date, Jimbo and 
Miwa [8], associated with an elementary auto-Backlund transformation 
for the continuous Chen-Lee-Liu hierarchy (cf. [21] )• Unfortunately, 
this system does not admit a complex/Hermitian conjugation reduction 
between u n and v n , so it cannot provide a proper discretization of the 
Chen-Lee-Liu equation. 
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• In the special case of a = — 5 and f3 = — 7 , (I2.13P reduces to a nontrivial 
lattice system 

' u n: t - ay ( 7 SI - u n v n ) (y 2 / + u n - ity) 1 (yu™ + (5rx„_i) 
i + b5 (yMn +1 + (5Mn) (<5 2 / + v n u n+ i) 1 ( 76 / - tyu n ) - (c - d)u n = 0 , 

< 

Vn,t - b5 (yd/ - v n u n ) (d 2 / + u„_iu n ) (dry + yv n _i) 

, + (8v n+ i + yty) (y 2 / + ma+i)' 1 (yd/ - w n u n ) + (c - d)v n = 0, 

while the Lax pair given by (12. 8 p and (12.111) becomes trivial in the 
naive limit a —> —6, f3 —> — y. This lattice system can be derived from 
a Backlund-Darboux transformation for the continuous NLS hierar¬ 
chy in the same manner as (12.131) is derived from a Backlund-Darboux 
transformation for the Ablowitz-Ladik hierarchy. By further setting 
y = d = ± 1 , we obtain a new proper space-discretization of the matrix 
NLS system as well as the matrix modified KdV system. We will dis¬ 
cuss it in a separate paper. 


Not all of the four parameters cy /3, y and d in (12.13p (or (12.140 1 are essential; 
the number of independent parameters can be reduced by applying a point 
transformation of the form: u n = /j,u n u' n , v n = v~ n v' n with nonzero constants 
fi and v. However, it then becomes difficult to consider limiting cases where 
one or more parameters tend to zero as mentioned above, so we prefer the 
“redundant” expression (12.13p (or (12.141) 1. 

By setting y = a, 5 = —j3 and a/3 = 1, (12.14p provides a new integrable 
space-discretization of the Chen-Lee-Liu system: 

a n,t 4“ a (Un ^n—l) (/ ^nMn— 1) (/ 3~ V n U n ^) 

I b (u n + 1 u n ) (/ T v n u n+ 1 ) (/ v n u n ) O, 

1 (2.15) 

Vn,t b (/ V n U n ) (/ T V n —i U n .) iyn ^n— l) 

+ a (/ + u n u n ) (/ - ry+iA,) -1 (u n+ i - v n ) = O. 

System (12.151) with b = —a* admits the Hermitian conjugation reduction 
v n = i v3 n ] in the case of a = b = i, we obtain a proper space-discretization 
of the Chen-Lee-Liu equation [Hj (see for the matrix case): 

+ (W„+1 - U n ) (/ + iujjlin+i) 1 (/ - i u\u n ) 

- (u n - u n _ i) (/ - 1 (/ + iu^M„) = o. 
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In the scalar case, this reads 


11 U r 


l^n,t 


1 + iM*M n+ i 

or equivalently, 


(^ra+1 U n ) 


1 + i |« ra | 2 

L - m*u n _i 


(u n u n — i) 0, 


iu n ,t T (^w+l “1“ U n —\ 2M n ) 


. (*4+1 - «n) < . K - *4-0 


Ml 


1 ~l“ hi n iin+i 


= 0 . 


lUnUn—l 


In addition, system (I2.15P with b = —a admits the matrix transpose reduction 
v n — u nC, where C is a constant skew-symmetric matrix |63j : in particular, 
in the case of a = 1 and b = — 1 , we obtain 


Un,t “h (Wn-i-i U n ) (d U n Cu n -\-l) (d U n CUn) 

+ (U n ~ U n - 1 ) (d - U^CUn- 1 ) 1 (d + U^Cu n ) = O, C 1 = ~C. 


In the scalar case, (12.13[) is a Hamiltonian system with an ultralocal (but 
noncanonical) Poisson structure. Indeed, it can be expressed as u n ^ = {u n , H} 
and v n) t = {u n , dd}, where the Hamiltonian and the Poisson brackets are given 
by 
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and 


{u m , U n } {M m , Vn\ 0, { U m , dmn (cV/3 U n V n ) ( 7 & U n V n ) , 

respectively. Here, S mn is the Kronecker delta, which should not be con¬ 
fused with the free parameter 6. This Hamiltonian structure encompasses 
the already known Hamiltonian structures for the simpler lattice systems in 
the scalar case (1^ - 12^1551 [53 1 [55]. It would be interesting to construct the 
corresponding classical r-matrix a la Sklyanin [ 651157 ]. 


3 Space-discrete massive Thirring model 

In this section, we describe how to discretize one of the two independent 
variables of the massive Thirring model in light-cone (or characteristic) co¬ 
ordinates [TTI 15o] . which is the first negative flow of the Chen-Lee-Liu hi¬ 
erarchy @31118] (see @6] 137] for the model in laboratory coordinates). With 
a slight abuse of terminology, we will call it a “space-discretization” of the 
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massive Thirring model, although it would be more natural and appropriate 
to use it for the massive Thirring model in laboratory coordinates. 

The generalized binary Backlund-Darboux transformation (12.6p with suit¬ 
ably defined intermediate potentials u m and v m can preserve the Ablowitz- 
Ladik spectral problem (12.11) form-invariant, as well as an infinite set of 
isospectral time-evolution equations such as (12.21) and (12.31) . Thus, the tem¬ 
poral Lax matrix M n given by (12.lip with (12.121) is not the only possible 
choice; it can be replaced with any M n -matrix corresponding to a flow of the 
Ablowitz-Ladik hierarchy, or equivalently, the continuous Chen-Lee-Liu hi¬ 
erarchy. To obtain a space-discrete analog of the massive Thirring model, we 
suppress the lattice index m for the Ablowitz-Ladik hierarchy and consider 
the isospectral time-evolution equation (12.31) written in the form: 
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Then, substituting the Lax pair given by (12.8p and 
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into the zero-curvature equation (12.101) . we obtain a set of six equations. A 
direct calculation shows that only four of them are independent, which can 
be presented as 


/5A (yQri- 1-1 C^Qn) T (a/3 7 $) ^nVn (flQn+l &Qn) O, 

(35 (aR n+ i - 7 R n ) - (a(3 - j5) v n - (5R n+1 - / 3R n ) u n v n = O, 

(a/3 y$) ^ n,t T ay ((3Qn- 1-1 5Q n / T u n (aR n+ i 7 -Rn) 

U n (O 3 //,; - \ (/n ■ 1 iSRnQn) O, 

(a/3 - yd) v n<t + ay (5R n+1 - (3R n ) + v n (y Q n+1 - aQ n ) v n 
k ( r )5Rn+lQn+l Ot(3RnQn) ^n O. 


By setting y = a = i and 8 = —/3 = —2i/A where A G R is a lattice param- 
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eter, (j3.ip reads 


(Qn+1 Qn) d" 2iu n U n V n (Qn+l d" Qn) Q? 

2 

(Rn+1 Rn) 2 iu n “h (Rn+1 d~ Rn) nVn Q? 


< 2 'U n ^ + i (Qn+1 d" Qn) 2 (Rn+1 Rn) ^n (3.2) 


^n (Rn+lQn+1 d" RnQn) Qj 


2^n,£ i (Rn+1 d" Rn) ^ (Qti+ 1 Qn) ^n 


d“ (Rn+lQn+1 d" RnQn) ^n Q* 


Note that the set of coefficients in (13.21) can be changed by rescaling the vari¬ 
ables. System (13.2p admits the Hermitian conjugation reduction R n = iQf, 
v n = iu^ , so we obtain a proper space-discretization of the massive Thirring 
model as 


( 2 


(Qn+1 Qn) d" 2lU n YLL n v) n (Qn+1 d“ Qn) O, 



(3.3) 


In the continuous space limit, (13. 3 li indeed reduces to the massive Thirring 
model in light-cone coordinates (see (3.48) in [44]): 

{ Q x + m — i uu'Q = O, 
u t + iQ — i uQ^Q = O. 

We remark that another lattice version of the massive Thirring model in 
light-cone coordinates was studied in [ 68 ], 

4 Concluding remarks 

In this paper, we have developed an effective approach for generating new in- 
tegrable lattice systems from Backlund-Darboux transformations for known 
integrable systems. The idea to interpret a Backlund-Darboux transforma¬ 
tion as a discrete spectral problem is already well-known (see, e.g., [SHE]); 
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the main new feature of our approach is to consider the intermediate po¬ 
tentials appearing in a binary Backlund-Darboux transformation as new de¬ 
pendent variables in the discrete spectral problem, up to a rescaling of the 
variables. As the name implies, the binary Backlund-Darboux transforma¬ 
tion is equivalent to the composition of two elementary Backlund-Darboux 
transformations in either order of the composition MUD]. In the example 
of the Ablowitz-Ladik hierarchy, u m and v m in (12,6j) appear in the decompo¬ 
sition of the binary Backlund-Darboux transformation into two permutablc 
elementary Backlund-Darboux transformations as |50j 

(Qm; R j m ) }■ (u m , *) t (Q m , Rm ) 


and 


(Qnm Rm) t (*, V m ) t (Q m , Rm ) 


up to a rescaling and a space translation, so they are indeed the intermediate 
potentials; a similar decomposition holds true for the continuous Chen-Lee- 
Lin hierarchy if one considers the pair of variables ( Q m , R m - 1 ) or (Q m -i, R m ) 
(cf. (12.2p or (I2.3p h Then, we associate the new discrete spectral problem with 
a suitable isospectral time-evolution equation on a case-by case consideration 
and obtain an evolutionary lattice system from the compatibility condition 
called the zero-curvature equation. We illustrated this approach by deriving 
the general lattice system (I2.14p . which involves the arbitrary parameters a, 
(3, 7 and 5 and includes the space-discrete Chen-Lee-Liu system (12.15p as 
a special case; the derivation of negative flows of the integrable hierarchy is 
easier and can be performed in a more systematic manner as is illustrated in 
section [3] This approach is quite useful for obtaining proper discretizations of 
integrable systems that admit the complex/Hermitian conjugation reduction 
between the two dependent variables, such as the NLS system, derivative NLS 
systems and their matrix generalizations; this is because the intermediate 
potentials generally inherit the internal symmetries of the original continuous 
system, which guarantees the feasibility of such a reduction. 

In the continuous case, by applying a nonlocal transformation of depen¬ 
dent variables 141114811711172 ] , we can transform the Chen-Lee-Liu equation to 
other derivative NLS equations such as the Kaup-Newell equation mi and 
the Ablowitz-Ramani-Segur (Gerdjikov-Ivanov) equation [ 73 ][ 72 ]; the same 
transformation applies to other flows of the integrable hierarchy including 
the massive Thirring model [33S3SEII52J (also see [56]). Note that the non¬ 
local quantity used in this transformation is given by the first (or second) 
component of the linear eigenfunction of the Lax-pair representation with 
the spectral parameter set equal to zero. In a similar manner, we can ob¬ 
tain new proper space-discretizations of other derivative NLS systems from 


14 







the space-discrete Chen-Lee-Liu system 02.151) by applying a discrete ana¬ 
log of the nonlocal transformation. Indeed, 02.71) with 02.8b|) implies that 
the quantity Z n := (— 5)~ n e^ lim ^ 0 C n ^i,n satisfies the following relation for 
any constant £: 

Z n -{-i (a(3I u n v n ) ( 7 <5/ -I- u n v n ) Z n . 


The time derivative of Z n for a suitably chosen £ can be obtained from 02.9|) 
with 02 . 111 ) and 02 . 121 ) as 

Z n>t = (a/3 - jS) [a /?7 (flyl - w n _iu n ) _1 + baS (aSI - Z n . 

Thus, the nonlocal transformation 


or its one-parameter generalization in the scalar case 

Qn ■ Zjf Un , f 'n ■ 7 ( Z n 


can be applied to 02. 14)1 : by setting 7 = a, 5 = —(3 and a/3 = 1, (12.14)) reduces 
to 02.15p and we obtain proper space-discretizations of other derivative NLS 
systems. 

In the original formulation of Backlund-Darboux transformations 
[3(71432] . the intermediate potentials such as u m and v m in 02 . 6 j) can be writ¬ 
ten explicitly in terms of the linear eigenfunctions at some fixed values of 
the spectral parameter (. Then, Backlund-Darboux transformations can 
be applied iteratively to obtain a sequence of new solutions of a nonlinear 
integrable system from its seed solution and the associated linear eigenfunc¬ 
tions; the final result does not depend on the order of applications. That is, 
Backlund-Darboux transformations with (generally) different values of the 
Backlund parameters are mutually commutative as long as boundary con¬ 
ditions are fixed appropriately; this fact can be understood intuitively by 
identifying Backlund-Darboux transformations as discrete-time flows that 
belong to the same integrable hierarchy. Without taking into account how 
to express the intermediate potentials in terms of the linear eigenfunctions 
of the Lax pair, the permutability condition for Backlund-Darboux transfor¬ 
mations results in a matrix re-factorization problem (see, e.g., [75HZ7]). In 
the case considered in this paper, it reads )Tl l7(H 174] 


L(U', V ; fi, v, £, rj)L(u, v; a, (3, 7, 6) = L(u' , v'\ a, (3, 7, 5)L(U, V; p, u, £, 77 ), 


(4.1) 
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where the matrix L(u,w,a, (3,^,8) is given by (j2.8p with the unnecessary 
lattice index n removed and the dependence on the spectral parameter ( 
suppressed. Considering the matrix inverse of L, we can rewrite (14.11) as US 

[ 73 ] 


L(u, v; 7, 8, a, fi)L(U', V'; f , 77 , /j, v) = L(U, V; f , 77 , //, i/)L(u', v'; 7, <5, a, /3). 


(4.2) 


Equation (14.11) or (j4.2j) can be naturally represented using an elementary 
quadrilateral and defines a parameter-dependent Yang-Baxter map 


(u,v,U',V) ^ (u!,v’]U,V), 


wherein the fields are assigned to the edges of the quadrilateral, instead of the 
vertices [771 - 180] . The Lax (or zero-curvature) representation for the Yang- 
Baxter map is given by (14.21) [76l 179] . It is not evident from the explicit form 
of the Lax matrix in (12.81) that (14.21) can be solved uniquely to provide the 
Yang-Baxter map, but the map itself can, in principle, be constructed using 
the definition of the intermediate potentials u and v in terms of the linear 
eigenfunctions. 
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